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The effect of a magnetic field on the dipole echo amplitude in glasses (at temperatures of about 
10 mK) induced by the dipole-dipole interaction of nuclear spins has been theoretically studied. It 
has been shown that a change in the positions of nuclear spins as a result of tunneling and their 
interaction with the external magnetic field Eh lead to a nonmonotonic magnetic field dependence 
of the dipole echo amplitude. The approximation that the nuclear dipole-dipole interaction energy 
Ed is much smaller than the Zeeman energy Eh has been found to be valid in the experimentally 
important cases. It has been shown that the dipole echo amplitude in this approximation may 
be described by a simple universal analytic function independent of the microscopic structure of 
the two-level systems. An excellent agreement of the theory with the experimental data has been 
obtained without fitting parameters (except for the unknown echo amplitude). 

PACS numbers: 61.43.Fs, 76.60.-k, 81.05.Kf 



I. INTRODUCTION 

Glasses at temperatures below 1 K are known to ex- 
hibit a number of universal properties almost indepen- 
dent of their composition and fundamentally different 
from the properties of similar crystals. These properties 
are traditionally described within the model of two-level 
systems (TLS's) [l|. One of these properties is a two- 
pulse electric dipole echo in glasses, which is the delayed 
response to two high-frequency electromagnetic pulses. 

When a glass is subjected to two rf pulses with a fre- 
quency of about 1 GHz separated by a time interval r 
much longer than the pulse length, one can observe a re- 
sponse in the glass polarization at time r after the second 
pulse. Especially interesting is the pronounced nonmono- 
tonic magnetic field dependence of the two-pulse echo 
amplitude observed at temperatures of 10 mK in the ab- 
sence of paramagnetic centers in the glass @ • The nature 
of this interesting phenomenon remained unclear until 
works |3j, U, Ifj, l(| , where this dependence was associated 
with the presence of atoms with a nuclear quadrupole 
moment in the glass. 

Owing to the interaction of the nuclear magnetic mo- 
ment with the external magnetic field (Zeeman interac- 
tion) and the nuclear quadrupole moment with the gra- 
dient of the internal electric field, the levels of the TLS 
split into two nearly identical series of levels. The char- 
acteristic energies of such a fine structure are 10 -21 erg 
(corresponding to a frequency of about 100 kHz), which 
are much lower than the typical energy of 50 mK of the 
TLS's specified by the 1 GHz frequency of the excita- 
tion pulses. Since the internal field gradient is different 
at different positions of the TLS's, the fine level struc- 
ture depends on the tunneling state of the TLS, which 
ultimately leads to the nonmonotonic magnetic field de- 
pendence of the dipole echo amplitude [!, H, Q . 



As a proof of this hypothesis, the results on measur- 
ing the two-pulse echo amplitude in glycerol (C3H8O3) 
were presented [J|. In this experiment the echo ampli- 
tude dependence on magnetic field increased by more 
than an order of magnitude under the replacement of 
hydrogen with zero nuclear quadrupole moment by deu- 
terium (whose nuclear spin is 1 so deuterium has a small 
quadrupole moment). 

The oscillations of the echo amplitude in a magnetic 
field were theoretically studied in [1, S B S Q • In par- 
ticular, the qualitative explanation of the observed effect 
was given in 0, Q . In 0, [|[ , numerical calculations 
of the dipole echo amplitude in glycerol in the magnetic 
field were performed and a good agreement with the ex- 
perimental data was obtained. 

In our previous work 0], we qualitatively and, in some 
cases, quantitatively compared our analytic expressions 
for the echo amplitude in glasses containing nuclear elec- 
tric quadrupole moments with the experimental data. 
However, the existence of small oscillations of the dipole 
echo amplitude in glasses containing only spherical nu- 
clei (without quadrupole moments), like in nondeuter- 
ated glycerol C3H8O3, remained unexplained until pa- 
pers [3, §• The presence of a small (0.037%) natural 
impurity of the 17 isotope, which also has a nuclear 
quadrupole moment, did not explain these oscillations. 

Fleishmann 0] suggested that the magnetic field de- 
pendence of the echo amplitude in the case of non- 
deuterated CsHgOs may be caused by the dipole-dipole 
interaction of the nuclear magnetic moments of hydrogen 
atoms (the spin and, respectively, the magnetic moment 
of the primary isotopes of carbon 12 C and oxygen 16 O are 
zero). This interaction (along with the Zeeman interac- 
tion) induces a fine structure of the levels of the TLS, 
which depends on the tunneling states of the systems. 

Bazrafshan et al. @ numerically calculated the echo 
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amplitude in deuterated glycerol C3D5H3O3 (taking into 
account all of the spins of hydrogen and the quadrupolc 
moments of deuterium) . It was assumed that the tunnel- 
ing of the two-level system is the rotation of the glycerol 
molecule as a whole. Although the theory was in a good 
agreement with the experiment in high fields and in semi- 
quantitative agreement in low magnetic fields (when the 
dipole-dipole interaction is responsible for the effect), the 
assumption of the rotational character of tunneling was 
not justified. 

This work is aimed at obtaining the analytical results 
for the magnetic field dependence of the dipole echo am- 
plitude in the absence of nonspherical nuclei in the glass, 
so that only the dipole-dipole interaction of the nuclear 
spins is taken into account, and at comparing the results 
with the experimental data on non-deuterated glycerol. 
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Here fii = //jJj/Jj is the operator of the magnetic mo- 
ment of the i-th nuclei, Jj is the operator of its spin, 

and H is the external magnetic field. The radius vector 
(1 2) 

Tij = r\j : from the i-th to the j-th nuclei may depend on 
the TLS position (1) or (2). The summation is performed 
over all tunneling nuclei. Consequently, 



Wj = 
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II. GENERAL THEORY OF THE MAGNETIC 
FIELD DEPENDENCE OF THE ECHO 
AMPLITUDE 

Parshin Q derived a general formula specifying the 
dipole echo amplitude in an arbitrary system, which con- 
sists of two identical sets of levels separated by an energy 
gap much larger than the splitting within the sets 
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Here, P oc ho is the dipole echo amplitude; V\ and V2 are 
the amplitudes of the first and second exciting electric 
pulses, respectively; E n are the fine-structure energy lev- 
els of the TLS, N is the total number of this levels, and 
ctnm is the matrix element for the transition between the 
n-th lower and m-th higher levels of the split TLS during 
the action of the excitation pulses. 

To use this formula, we have to find the fine-structure 
energy levels and the matrix elements of the transition 
in an arbitrary magnetic field. For that, let us write the 
Hamiltonian of the system including, first, the tunneling 
of the TLS atoms and, second, the energies of the nuclear 
spins, which form the fine level structure: 



H tls ®lj+% ls ®Wj+a zMs ®Vj + {d as F)®ij. (2) 



Here Ht\ s is the Hamiltonian of the TLS in the coordinate 
representation (excluding the fine structure); Lis and lj 
are the identity matrices in the spaces of the TLS (2 x 2) 
and the nuclear spin projections, respectively; cr z ,ti s is 
the Pauli matrix for the two-level system; dti s = d<j Z) ti s 
is the operator of the electric dipole moment d of the 
TLS; and F is the external (excitation) electric field. 

The spin operators Wj and Vj are, respectively, the 
symmetric and antisymmetric (with respect to the TLS 



Note that the first term of spin Hamiltonian ([3]) , which 
is associated with the external magnetic field (the Zee- 
man part), does not include vectors that change under 
tunneling of the TLS and, therefore, does not contribute 
to the antisymmetric Vj part of the Hamiltonian. The 
dipole-dipole interaction (the second term of Hamiltonian 
([3])) enters both operators Wj and Vj. 

To find the quantities associated with the fine level 
structure, let us rewrite Hamiltonian @ in the basis con- 
sisting of the stationary wavefunctions of the bottom and 
top levels of the TLS and the eigenfunctions of the oper- 
ator Wj for the spin variables: 
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Here, A is the difference between the minima of the 
double-well potential, A is the tunneling amplitude of 
the initial (unsplit) TLS, and E = ^A 2 + A 2 is the to- 
tal energy of the TLS (without the inclusion of the fine- 
structure effects). The tilded operators W j and V j cor- 
respond to the choice of the eigenfunctions of Wj as the 
basis. In this basis, Wj is a diagonal matrix. 

We assume that the characteristic tunneling lengths of 
the nuclei are much shorter than the characteristic inter- 
atomic distances and, consequently, V j is much smaller 

than W j and may be taken into account perturbativcly. 

Retaining only the terms proportional to V j in the tran- 
sition matrix elements, we obtain the dipole echo ampli- 
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matrix form of Wj is 
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Here, e nm = [Wj) - ( Wj) = E n - E m is the 

\ / nn V / mm 

difference between the fine-structure levels of the TLS. 

Expression ([6]) allows one to calculate the echo am- 
plitude from one TLS. To find the observed dipole echo 
amplitude of the entire sample, it is necessary to perform 
the summation over all two-level systems. 

To obtain the expression for the distance between the 
fine-structure levels e nm , we have to diagonalize the ma- 
trix Wj, which cannot be done analytically even in the 
simplest case of two interacting spins 1/2 (a 4x4 matrix). 
Thus, to go further, we need to perform some simplifica- 
tions. 

The nonmonotonic magnetic field dependence of the 
dipole echo amplitude in C3H8O3 was experimentally 
observed in the magnetic fields H ~ 50 G, which cor- 
responds to the Zeeman energy Eh — 0.7 x 10 -21 erg. 
The energy of the dipole-dipole interaction between two 
hydrogen nuclei at a distance of 1.9 A (the minimum dis- 
tance between hydrogen atoms in glycerol) may be esti- 
mated as Ed — 1.5 x 10 -22 erg, which is almost an order 
of magnitude smaller than the Zeeman energy. 



III. WEAK DIPOLE-DIPOLE INTERACTION 
APPROXIMATION 

In the Ed <C Eh approximation, we neglect the dipole- 
dipole interaction in the expression for Wj. However, 
this interaction should be retained in Vj, which does not 
contain terms of the order of Eh- 

In this approximation, the eigenfunctions of Wj are 
the states with definite projections of all nuclear spins on 
the direction of the magnetic field. In this case, only the 
matrix elements (Vj) nm of the transition with a change 
of no more then two spin projections are nonzero. This 
allows us to reduce the problem of an arbitrary number of 
dipole-dipole interacting nuclear spins to the problem of 
one pair of spins. We consider only the pairs of identical 
spin 1/2 nuclei. The operator Wj may be then written 
in the form 



Wj = (j,H(3 1 +3 2 ) 
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where /i is the nuclear magnetic moment and Ji and 
J2 are the spin operators of the interacting nuclei. The 
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It is easy to see that, in this case, there are four energy 
differences e nm equal to 2fiH and one difference £14 = 
4/j.H (£23 is zero and does not contribute to the echo 
amplitude) . 

To express the echo amplitude explicitly with the use 
of Eq. ©, we need to calculate also the average squares 
of the matrix elements |(Vj) nm | 2 . Taking into account 

the smallness of the tunneling length, the V j operator in 
the case of two nuclei can be represented as 



Vj = Sr 12 
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Expressing the gradient explicitly and taking out Sr 12, 
ri2 and ^1,2 as the common factors, we may rewrite Eq. 
© in the form 
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Here, the Greek subscripts numerate the Cartesian co- 
ordinates; ai <a and a2,p are the Pauli matrices for the 
first and second spin, respectively; iV a ^(n r , n$ r ) is the 
symmetric second-rank tensor, which depends merely 
on the directions of the vectors r (n r = r/r) and 6r 
(ng r = Sr/Sr). 

To find the sought expression for |(Vj)„ m | 2 , one has to 
express one matrix element from Eq. (|10[) and average its 
absolute value square over the directions of the vectors n r 
and n$ r (it is easy to see that the averaging affects only 
the tensor N a p). For this purpose, we introduce a new 

numeration of the elements of V j : ( Vj) 
Here, the subscripts p,q = ±1 and r, s 
the states of the first and second spins, respectively. Cor- 
respondingly, n = (pr) and m = {qs). In this notation, 

§ r 2 

\{Vj)( pr )( qs )\ 2 = M 4 — X 

(N a pNs C ). (11) 

Here, ci^p? is the pq element of the Pauli matrix a a 
acting on the wavefunctions of the first spin. tJ2,i3,rs is 
the same for the second spin. 

The isotropy of the glass allows us to express the 
symmetric tensor (N a pN^) in terms of two constants, 
(NapNtf) = ASapS^ + B(S ai S K + SacSpt). However, 
it can be shown by direct calculation that the quantity 

CT 2,/3,rs^i,£, M ^I,C ) rs ( W^£C = \ a i,vq ' a 2,rs\ 2 , which 
corresponds to the constant A, includes only the diagonal 
elements and the elements corresponding to the transi- 
tion between the degenerate second and third levels of 
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FIG. 1: Experimental data on glycerol 0,[T(J at (a) r = 3.5 fis 
and (b) 8 /is in comparison with Eq. (|12[1 with fi taken equal 
to the proton magnetic moment. 



Wj. Neither of them contributes to the echo amplitude 

and, consequently, all significant matrix elements of V j 
may be expressed in terms of a single constant B. 

Performing the calculations, we arrive at the final ex- 
pression for the dipole echo amplitude in the form 
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where C is a constant independent of the magnetic field. 

Note that expression (fT2")l was derived without any as- 
sumptions about the relation between n r and n,5 r and, 
consequently, is absolutely independent of the tunneling 
geometry of the atoms, that is, of the microscopic struc- 
ture of the TLS. 

Let us compare Eq. (fl"2)) with the measurements [1, [l(| 
of the magnetic field dependence of the dipole echo am- 
plitude in glycerol (see figure). Clearly, the analytical 
curve with [i equal to the magnetic moment of a proton 
describes the experimental data sufficiently well without 



any fitting parameters (except for the constant C, which 
specifies the vertical scale) until the region of small mag- 
netic fields (< lmT), where the conditions Eh <C Ed 
seemingly breaks. 



A. Applicability Conditions 

The nonmonotonic magnetic field dependence of the 
dipole echo amplitude is determined by the term 
sin 4 (e nm r/2fi) in Eq. ©. This corresponds to the char- 
acteristic energies e nm ~ 2K/t. In the case of Ed <C 2H/r 
(or, equivalently, r <C 2h/E ( i), we may expect that the 
condition Ed <C Eh holds for the dominant part of the 
dependence. Otherwise, the exact inclusion of the dipolc- 
dipolc interaction is required; in this case, in particular, 
the problem cannot be reduced to the problem of two tun- 
neling nuclei. This agrees well with the fact that the the- 
ory better describes the experimental data for r = 3.5 fis 
than for r = 8 /j,s (see figure) in the low-field region. 

Note that the horizontal scale of the magnetic field 
dependence of the echo amplitude in the case under con- 
sideration is determined by the time t. Thus, the curve 
Pccho(H) should scale in the horizontal axis as 1/r with 
a change in the time interval between the pulses. In the 
opposite case of r > 2H/Ed, the horizontal scale of the 
dependence is determined by the relation Ed — Eh and 
is independent of the time interval between the pulses. 
This may be used as an experimental criterion of the 
smallness of the dipole-dipole interaction. 

To conclude, in the experimentally observed case of 
r <C 2h/Ed, the magnetic field dependence of the dipole 
echo amplitude induced by the dipole-dipole interaction 
of the nuclei of the same kind has a universal character 
independent of the microscopic structure of the TLS. It 
is well described by a simple analytic expression. A good 
agreement of the theory and experiment on the dipole 
echo in nondeuterated glycerol was obtained without any 
fitting parameter (except for the echo amplitude scale). 
At larger values of the time interval between the excita- 
tion pulses, t > 2h/Ed, the universal form of the depen- 
dence breaks, which may allow one to obtain information 
on the microscopic structure of the two-level systems with 
the use of numerical analysis [lll |. 
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